arXiv:1505.07685vl [cond-mat.str-el] 28 May 2015 


Incommensurate spin density wave at a ferromagnetic quantum critical point in a 

three-dimensional parabolic semimetal 

James M. Murray/ Oskar Vafek/ and Leon Balents^ 

^National High Magnetic Field Lab and Department of Physics, 

Florida State University, Tallahassee 32310, USA 
^Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 

(Dated: May 29, 2015) 

We explore the ferromagnetic quantum critical point in a three-dimensional semimetallic system 
with upward- and downward-dispersing bands touching at the Fermi level. Evaluating the static 
spin susceptibility to leading order in the coupling between the fermions and the fluctuating ferro¬ 
magnetic order parameter, we find that the ferromagnetic quantum critical point is masked by an 
incommensurate, longitudinal spin density wave phase. We first analyze an idealized model which, 
despite having strong spin-orbit coupling, still possesses 0(3) rotational symmetry generated by the 
total angular momentum operator. In this case, the direction of the incommensurate spin density 
wave propagation can point anywhere, while the magnetic moment is aligned along the direction of 
propagation. Including symmetry-allowed anisotropies in the fermion dispersion and the coupling 
to the order parameter field, however, the ordering wavevector instead breaks a discrete symmetry 
and aligns along either the [111] or [100] direction, depending on the signs and magnitudes of these 
two types of anisotropy. 

PACS numbers: 


For decades it has been understood that electron in¬ 
teractions may lead to interesting consequences at low 
energy scales in three-dimensional (3D) semimetallic sys¬ 
tems, in which spin-orbit coupling combines with cu¬ 
bic crystalline symmetries, leading to a band degener¬ 
acy at the Fermi energy.^HSl Such a band structure has 
long been known to occur in HgTe and a-Sn (gray tin), 
and has possibly been discovered more recently in the py- 
rochlore Pr 2 lr 207 .^ 2 Hni Since the low-energy effective the¬ 
ory for these systems was first developed by Luttingeili^, 
it has been argued using various approaches that in¬ 
teractions c ould lead to an excitonic instability at low 
temperature^^, or that the system may be described 
by an exotic non-Fermi liquid phase, characterized by 
nontrivial power-law scaling of various physical quanti¬ 
ties with temperature.f^lf^ In more recent work, the prop¬ 
erties of such a system near an Ising antiferromagnetic 
quantum critical point were explored, and the critical 
theory was found to be governed by unusual critical ex¬ 
ponents and emergent spatial anisotropy of the fermion 
dispersion.l^ The theory describing the quantum phase 
transition into an insulating nematic phase has also been 
developed recently.l^ 

In this work we explore the fate of a 3D semimetal 
in the vicinity of a ferromagnetic (FM) quantum criti¬ 
cal point. Working at zero temperature and approaching 
the quantum critical point from the paramagnetic side, 
we find that it is unstable toward an incommensurate 
spin density wave (SDW) phase. In the most symmetric 
version of the theory, there is a combined 0(3) rotational 
symmetry for spin and spatial degrees of freedom, which 
is spontaneously broken by the SDW wavevector. The 
0(3) symmetry is reduced to a discrete symmetry by ei¬ 
ther of two sources of anisotropy: the Yukawa coupling 
of the fermions to the fluctuating FM field, or the disper¬ 


sion of the fermions themselves. (Throughout this work, 
by anisotropy we shall mean terms preserving the lattice 
point group symmetries, but not 0(3) symmetry.) Ei¬ 
ther of these sources of anisotropy reduces the 0(3) sym¬ 
metry breaking to a discrete symmetry breaking, with 
the ordering wavevector lying along one of several high- 
symmetry directions. In contrast to the more typical 
case of an incommensurate SDW driven by Fermi surface 
nesting, the incommensurability of the SDW found in 
our work is unrelated to Fermi surface or doping effects. 
Rather, the ordering wavevector depends on nonuniver- 
sal parameters such as the temperature and strength of 
the boson-fermion coupling. 

In the following section we introduce the model Hamil¬ 
tonian for a 3D parabolic semimetal coupled to a fluctu¬ 
ating FM order parameter field and discuss its relevant 
symmetries. In Section we use this model to study 
the effects of fermionic fluctuations on the bosonic order 
parameter field, beginning with the 0(3)-symmetric case 
and then moving on to the anisotropic case, with details 
of these calculations provided in the Appendix. In both 
cases we obtain a negative contribution to the boson self¬ 
energy that is linear in momentum, leading to an instabil¬ 
ity to an incommensurate SDW phase. In Section [ITT] we 
investigate the possibility of a fluctuation-induced first- 
order phase transition directly from the paramagnetic to 
the FM phase, evading both the FM quantum critical 
point and the incommensurate SDW phase. Comput¬ 
ing the free energy of the FM order parameter, we find 
that the sign of the fluctuation-induced nonanalytic term 
is such that no first-order transition results, in contrast 
with the metallic case in which there is a nonvanishing 
density of fermions at the Fermi energy.l^sHIZl Finally, in 
Section |IV] we discuss the results and their possible rele¬ 
vance to real materials. 
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The five T matrices satisfy {rijTj} = 25ij, and one can 
build ten additional matrices using Tij = To¬ 

gether with the unit matrix these make a complete ba¬ 
sis of 4 X 4 matrices. We use the T-matrix representa¬ 
tion, spin matrices Si, and basis functions (ii(k) given 
in Ref. UHl For simplicity we assume particle-hole sym¬ 
metry (cq = 0), which was found to emerge in the low- 
energy theory from the renormalization group calcula¬ 
tions of Refs. I5I7I 

The bare fermion Green function is 


FIG. 1: Schematic phase diagram showing the onset of incom¬ 
mensurate spin density wave (ISDW) order that preempts the k) = 

quantum phase transition to a ferromagnetic (FM) phase at ~ 

r < 0. The nature of the phase transition between the ISDW 
and FM phases (not shown) is not specified by our theory. 

where 


-inU-Hoik) 


( 4 ) 


I. THE MODEL 

The following action describes a theory of fermions 
coupled to a fluctuating FM order parameter (j) at tem¬ 
perature T — 0: 


S = J dr +'Ho{-i'V)]ilj + (f ■ ^-0 


■0^ (Ml cos a + M 2 sin Oi)tp • 0] >. 


( 1 ) 


y]v 


The fermion fields tp have 4iV components, where = I 
in the physical case. The action ([^ is similar to the one 
studied in the Ising antiferromagnetic cas^, but with 
some important differences. Most obviously, the FM or¬ 
der parameter 0 is a vector, rather than a scalar as in 
the AF case. The mass tq = (0, 0)])^^ of this field 

is tuned to zero at the FM quantum critical point , as 
shown in Figure Finally, as pointed out previousljD^, 
there are two symmetry-allowed terms through which FM 
order can couple to the fermions. Hence u describes the 
overall strength of the coupling between the fermions and 
the field 0, while a parametrizes the relative strength of 
the two allowed couplings. The 4x4 matrices ap pear ing 
in Q are related to the S' = 3/2 spin matricei^^ by 
Ml = (S„ Sy, S,) and M 2 = {Si S^, S^). 

The noninteracting part of the fermion Hamiltonian 
appearing in Q is the same as that appearing in Refs. 0 


'Hoik) = cok^ + Cl X! + C2 X! (2) 


n—1 


n—4 


where 


(ii(k) = Viikyk^, d2{k) = dspk) = VSk^ky, 

Md \ (3) 

d4(k) = \{kl- kl), d,{k) = -{2kl -kl- kl). 


Ek = \Jcl [dl (k) -b dl (k) -b dl (k)] -b (k) -b d^ (k)]. (5) 


The two-fold degenerate energy bands thus have disper¬ 
sion ±Ek, with a quadratic band touching point at k = 0. 
In the isotropic case in which ci = C 2 , it is straightfor¬ 
ward to show that [L -|- S, 'Ho(k)] = 0, where L = r x k, 
[ri,kj] = iSij, and L -|- S is the generator of rotations. 
In this case the fermionic theory has complete rotational 
invariance, and the dispersion from (1^ becomes simply 
Ek = cik^. 


II. FERROMAGNETIC POLARIZATION 
TENSOR AND SDW INSTABILITY 


In this section we calculate the bosonic polarization 
tensor, which describes the self-energy of the bosonic field 
due to damping of spin fluctuations by fermionic excita¬ 
tions. We begin with the isotropic case and then consider 
the anisotropic case in the following subsections. 


A. Polarization tensor: isotropic case 


Let us evaluate the polarization function, first in the 
relatively simple case where c; = C 2 . Subtracting off the 
UV-divergent contribution n^ (0,0) = —Siju‘^Ac/{2Tr'^ci) 
(where is the ultraviolet momentum cutoff), which can 
be absorbed into a redefinition of the bosonic mass term 
by letting r = ro -t- H^ (0, 0), one obtains the following: 
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n;^(*f],q)-n*;(o,o) = 


2 poc> 


dPk 


dw 


Tr[G'o(iw + in/2, k + q/2)M*(a)Go(ia; — in/2, k — q/2)M'^ (a)] 




N J (27r)3 2tt 

- Tr[Go(iw, k)M\a)Go{iuj, k)M^'(a)] 
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fi (fl) Tr [M\a)M^{a)] +5^/3 {§) (q • A"* • q)(q • A" • q) + /4 (^) q • A^A" 


( 6 ) 


q 


+ /5(^)Tr[A" 




Tr [r„M*(a)r„M^'(a)] , 


where we have set ci = C 2 = 1, and defined q = q/|q|, 
M^{a) = {Ml cos a + M^sina), and q • A"* • q = 
The matrices A"* (not to be confused 
with the momentum cutoff A^) are proportional to the 
five symmetric Gell-Mann matrices and are given in the 
Appendix, where we also provide the scaling functions 
fi{in/q^). If it is further assumed that a = 0, so that 
the magnetic coupling has complete 0(3) symmetry, then 
([^ becomes the following (noting that the traces in ([^ 
each give a factor of N): 

n;^(*o,q)-n;^(o,o)^^ 



u |q| 


6/3 (f) + y/4 (f) 


,,,, , 15 
„2 ) + 2 


/s ( 


qiQ] 


Including the bare kinetic terms, 
stiffness is thus 


the dressed bosonic 


^d(7) FIG. 2: The scaling functions appearing in the dressed prop¬ 
agator for the FM field, given by The solid line shows 
Fi (ifl/q^), and indicates a finite energy cost for transverse 
spin fluctuations at all frequencies. The dashed line shows 
F 2 (ifl/q^), which describes the energy cost of longitudinal 
spin fluctuations and becomes negative for |f2|/q^ < 0.30, 
indicating a phase instability. 


[G4tn,q)]-; = [G6')(*II,q)] 




+ M^|q| 


Fi (f) {S^J - mj) + F, (f) 


1 ( 8 ) 


where we have introduced 

Fi=5/1 (f) - h (i) - h (f)+f /5 {§) , 

F 2 = 5/1 (fl) + 5/3 (fl) + f/4 {§) + fh (f) . ^ 


The factor in (|^ proportional to {5ij — qiqj) describes 
the dynamics of transverse spin fluctuations, while the 
factor proportional to qiqj describes longitudinal fluctu¬ 
ations. Figure|^shows the scaling functions appearing in 
([^ . While the transverse part remains positive for all val¬ 
ues of |n|/q^, implying that these spin fluctuations have 
a nonzero energy cost at all frequencies, the longitudinal 
part becomes negative at frequencies |n| < 0.42q^, im¬ 
plying a vanishing energy cost for longitudinal spin fluc¬ 
tuations at sufficiently small frequencies and momenta. 

To O(q^), the static part of bare boson propagator has 
the following analytic form: 


[G^°^(0,q)]i/ = r6,j + q^[u^(d,j - qiqj) + v'^q.qj], (10) 


with stability of the bare theory requiring that uj > 0 
and > 0 for r > 0. From ([^ it is evident that if 
one fixes 11 = 0 and decreases momentum |q| with r 
sufficiently small, there will be a critical value of |q| where 
the negative correction M^F 2 ( 0 )|q| will overtake the bare 
term. The consequences of this become evident upon 
inverting the expression ([^ at 12 = 0 to obtain the static 
spin susceptibility: 


Xii(0,q) 


4 . + yq2 + y(o)„2|q|(^b Fdj) 
1 

r + w|q2 -I- F'2(0)u^|q| ’ 


( 11 ) 


where the constants are F’i(O) « 0.0362, and F2(0) « 
—0.0042. Equation 0 makes it clear that the vanish¬ 
ing longitudinal term for sufficiently small r corresponds 
to an instability toward a spin density wave phase at 
finite wavevector IQI = \F2{0)\u'^/2v2- This incommen¬ 
surate instability is similar to the “spiral SDW” instabil¬ 
ity that has been noted in previous theories of itinerant 
fermion systems near a quantum critical point described 
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using the spin fermion modeP^I^, and to the behavior 
recently observed in the metallic ferromagnet PrPtAlf^ 
In the isotropic case that we have so far considered, the 
wavevector Q of the SDW order could be along any di¬ 
rection, and thus breaks an 0(3) symmetry, with the 
spins aligning along (since the SDW is longitudinal) the 
direction of Q. The fact that the spins cannot inde¬ 
pendently choose a direction along which to align is a 
direct consequence of the spin-orbit coupling built into 
this model, and distinguishes this theory from the spin 
fermion theory of itinerant ferromagnets, for which the 
SU(2) spin rotation symmetry is independent of spatial 
rotations. Finally, we note that while our calculations 
are performed at T = 0, one can see that at T > 0 the 
temperature dependence will enter through the scaling 
function F 2 . This leads to temperature dependence of 
the ordering wavevector, which is an unusual feature in 
an SDW system. 


B. Polarization tensor: Yukawa anisotropy 


Before proceeding to the fully anisotropic case with 
both Cl 7 ^ C 2 and a 7 ^ 0 , in this section we investigate the 
case in which the Yukawa coupling is anisotropic {a 7 ^ 0) 
but the fermion dispersion remains isotropic (ci = 02 ). 
In this case the summations in may be performed to 
obtain a polarization tensor of the following form: 


n;^(*D,q)-n;^(o,o) = 


u^lql 


hi{'^,a) + h2C^,a)qi +h3(p,a)g. 


*^12) 


h4C3,a) + h5{^^,a){qi + q.) 




The scaling functions hi are linear combinations of the 
fi appearing in §, with coefficients depending on the 


Yukawa angle a. 

In order to determine the nature of the SDW instability 
in the anisotropic case, we can obtain the eigenvalues and 
eigenvectors of the polarization tensor for these two cases 
from the following equation: 


I 

M^lql 


n0(o,q) -110(0,0) ?/’0(q) 


= 7r0(q)i/:0(q). (13) 


For q ^ (1,0,0), the eigenvectors are aligned with 
the principal axes. In this case the longitudinal mode, 
for which = ( 1 , 0 , 0 ), has eigenvalue 7 r^^^(q) = 

hi + h 2 + h^, while the two degenerate transverse modes 
have 7 r^^^(q) = hi. For q ^ (1,1,1), on the other 
hand, the longitudinal mode has eigenvalue 7 r^^^(q) = 
hi + |/i 2 + |/i 3 -f h 4 -I- |h 5 , while the transverse modes 

have eigenvalues 7 r^^^(q) — hi + 7^h2 + ^h^. By tracking 
the minimum eigenvalue as a function of q, as shown in 
Figurej^ we find that the ordering occurs along [111] and 
is longitudinal for all values of the Yukawa parameter in 
the range 0 < a < n/2. 


C. Polarization tensor: fully anisotropic case 


Considering now the more general case in which both 
Cl 7 ^ C 2 , and a 7 ^ 0 , in this subsection we show the way in 
which these two types of anisotropy favor ordering along 
certain crystalline axes. The polarization tensor for the 
FM field (p in this general case is given by 


Bj 


(*D,q) -n;^(0,0) = (f ,q) Tr [M\a)M^(a)] (f ,q) Tr [M\a)T^M^ (a)] 


(14) 


+ E Tr [M\a)M^a)r^] + ^ (p,q) Tr (a)] [. 


The scaling functions appearing in (14) are given in 
the Appendix, Equations (A2)-(A4). Introducing the 
anisotropy parameter 6 for the fermion dispersion by let¬ 
ting ci ^2 = 1 ± these scaling functions can then be 
obtained perturbatively in 6, as shown in the Appendix. 


As pointed out in the two previous subsections, the 
behavior of spin excitations becomes increasingly domi¬ 
nated by the contribution from n 0 (iD,q) at small mo¬ 


menta and frequencies. Due to the anisotropy, cer¬ 
tain ordering wavevectors will be favored when either of 
these parameters is nonzero. In order to determine the 
wavevector of the SDW instability, we once again inves¬ 
tigate the minimum eigenvalue 7 r 0 (q) from (13). This 
quantity is plotted in Figure]^ as a function of the angle 
of q. We find that the ordering occurs along [111] when 
(5 > 0, while (5 < 0 favors ordering along [100]. With such 
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FIG. 3: (a) The polarization tensor eigenvalues as a function 
of the Yukawa angle a with 5 = 0. Blue (red) lines correspond 
to longitudinal (transverse) modes, while solid (dashed) lines 
have q ~ [100] ([111]). The lowest eigenvalue occurs for the 
longitudinal mode with q ~ [HI]- (b)-(d) The minimum 

eigenvalues 7r,^(q). The direction of the ordering wavevector 
for the incommensurate SDW will correspond to one of the 
potential minima, shown as red points, (b) Anisotropy in the 
Yukawa coupling leads to minima along [111] and equivalent 
directions, (c) The positive fermion dispersion anisotropy pa¬ 
rameter 5 > 0 leads to minima along [111] and equivalent 
directions, (d) 5 < 0 leads to minima along [100] and equiva¬ 
lent directions. 


anisotropy taken into account, the ordered state itself is 
now an incommensurate SDW with wavevector whose di¬ 
rection is locked to the crystalline axes, and both the zero 
temperature and T > 0 transitions will no longer be in 
the 0(3) ferromagnetic class as in the naive unrenormal¬ 
ized theory, but must be reconsidered in this light. For 
brevity and to maintain focus we do not do this here. 

One can also ask whether the incommensurate SDW 
instability survives in the case that 5 is not perturbatively 
small. In this case the polarization tensor can be evalu¬ 
ated by directly integrating the scaling functions appear¬ 
ing in ( [I^ numerically at D = 0 for a given value of 5. As 
shown in Figure]^ the instability occurs for most values 
of 5 when the Yukawa parameter lies in the physically ex¬ 
pected range 0 < a < 7r/2. The exception occurs in the 
case where 5 approaches 1, or equivalently for ci ^ C 2 . 
The physical case, however, is more likely in the oppo¬ 
site limit of 5 < 0 (or equivalently ci < C 2 ), which is the 




FIG. 4: Minimum eigenvalue of the polarization tensor as a 
function of the dispersion anisotropy, with the wavevector q 
along the [100] direction (blue) and [111] direction (orange). 
The incommensurate SDW instability occurs at the wavevec¬ 
tor for which the eigenvalue 7r,^(q) is most negative. The 
angle parametrizing the Yukawa interaction is set to a = 0 
and Q = 7r/3 in (a) and (b), respectively. 


case appearing most generically f or sim ple tight-binding 
models on the pyrochlore lattice.l^^EH Accordingly, the 
case Cl ^ C 2 was considered in the anisotropic RG the¬ 
ory of Ref. m Also note that we have calculated 7r,^(q) 
only for wavevector directions q ~ [100] and q ^ [HI]- 
Based on the results for small 5 shown in Figure and 
general symmetry considerations, it is natural to expect 
that the ordering will occur along one of these directions, 
though we have not shown this rigorously. 


III. FREE ENERGY EXPANSION 

Previous investigations of itinerant fermionic systems 
near an FM quantum critical point have found that the 
quantum phase transition can be precluded by either of 
two types of instability.!^ The first, in which fluctuations 
contribute a negative term to the bosonic self-energy, 
leading to an SDW instability at finite wavevector ]Q], 
was already discussed in the previous section. The sec¬ 
ond is a fluctuation-induced first-order phase transition 
into the FM phase, the possibility of which can be in¬ 
ferred by the presence of a nonanalytic term in the order 
parameter in the free energy expansion.^^!!!^ In this sec¬ 
tion we derive the Landau free energy for the FM order 
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parameter, showing that, although a fluctuation-induced 
nonanalytic term is present, its sign is such that no first- 
order instability results. (The absence of such a transi¬ 
tion has also been found in the Ising antiferromagnetic 
case.l^ 

The following effective action is obtained by integrat¬ 
ing out the fermions from the original action 0 : 


^eff -2 


dr / 


[G, 


(O)i-l 




— Tr In 


Go- - ^M(a) 


(15) 


The trace in the second term is a trace over matrix el¬ 
ements, as well as over frequency and momentum. The 
quartic term ~ 5|(()|^, with 6 > 0, is included since it is 
allowed by symmetry and arises in the low-energy effec¬ 
tive theory from integrating out high-energy modes. Due 
to the gapless fermionic excitations, one cannot simply 
follow HertJ^ and expand the logarithm in powers of 
(j) to obtain an effective theory for the order parameter 
field. Rather, a useful hrst step is to differentiate (15) 
with respect to the order parameter. And because we 
are interested in uniform states, we can evaluate the re¬ 
sult at zero external momentum and frequency, so one 
has = rSij, i.e. the bare mass of the bosonic 

field. The result is 


SF[$] SS,s 


Scj^i S(j)^ 




y2 , U f dil 




d^q 

(^ 


(16) 


X Tr 


Go-^(*fl,q)-^M(a).</) 


-1 


M\a) 


If we take the matrix inverse and specialize to the 
isotropic case in which a = 0, this becomes 


SF[(j)] V? 


dVt 

27r 


d^q 

(27r)3' 


u 




u 

~n' 


/ r2.q 


bcl?(j>i 

f - /(*^^,q,o) 


(17) 


where f{i^, q, (j)) is a complicated function, and in the 
second line we have added and subtracted the UV- 
divergent piece. By rescaling fl —>■ q —>■ 

in the second integral, and taking advantage of the fact 
that—due to rotational symmetry—the integral must 
only depend on the magnitude of (j) but not its direction, 
(ITtI) becomes 


sm 

S(j)i 


= r + b(f^ + - f{in, q, 0) 



(18) 


[/ (*^^,q, 1) - /(*^^,q,o)], 


from which the free energy is 


F[$] = a|(^p -b c|(^|G2 + 


(19) 


By evaluating the integral in (18) numerically, the second 
constant appearing in (19) is c « 0.019u^/A^. Since c > 
0, we conclude that the nonanalytic term does not lead 
to a first-order phase transition into the FM phase, and 
that the FM quantum critical point is preempted by an 
incommensurate SDW via the mechanism discussed in 
Section im 


IV. DISCUSSION 

In this work we have shown that the quantum critical 
point in a 3D parabolic semimetal near the onset of FM 
order is preempted by an instability to an incommensu¬ 
rate, longitudinal SDW phase. While the wavevector Q 
characterizing the SDW order breaks a continuous sym¬ 
metry in the idealized, 0(3)-symmetric version of the 
theory, including anisotropy terms allowed by crystalline 
symmetry leads to a discrete symmetry breaking, with 
Q along either the [111] or [100] direction, depending on 
the values of the anisotropy parameters as shown in Fig¬ 
ure [^ We have shown that the instability is present over 
a broad range of anisotropies, apart from the extremely 
anisotropic case where Ci 3> C 2 . Investigating this limit 
and its RG flows further, along the lines of Ref. [T] could 
be an interesting direction for future work. Finally, we 
have also investigated the possibility of a fluctuation- 
induced first-order transition into the FM phase, find¬ 
ing that, although a nonanalytic term in the free energy 
does exist at zero temperature, its sign is such that no 
first-order instability results. 

While in this work we have not focussed on the role 
played by Coulomb interaction, some general statements 
about its effects can be made. As pointed out in previ¬ 
ous wor kP3, Coulomb interaction tends to decrease the 
effective mass difference between the electron and hole 
bands, as well as the rotational anisotropy in the fermion 
dispersion. Because the SDW instability dominates the 
low-energy physics, the effects of Coulomb interactions 
will be secondary near the FM quantum critical point. 
(A possible exception is in the limit ci 3> C 2 , where— 
as discussed above—the SDW instability is absent. In 
this case we cannot rule out the possibility that the com¬ 
bined effects of Coulomb interactions and spin fluctua¬ 
tions could in principle lead to a stable anisotropic hxed 
point such as the one found in Ref. 0 ) 

While there is no obvious existing material described 
by the critical theory developed in this paper, it is possi¬ 
ble that the vicinity of a (preempted) FM quantum crit¬ 
ical point could be reached by doping existing parabolic 
semimetals such as HgTe or a-Sn with magnetic impu¬ 
rities such as Mn. Half-Heusler and Heusler materi als 
can exhibit narrow-gap or semimetallic behavioil^Sl^ as 
well as ferromagnetic critical point^^, although to our 
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knowledge the combination of these two conditions has 
not been shown to exist in any one material. 
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Appendix A 


In this Appendix we derive the scaling functions determining the polarization tensor given in the main text. In 
order to compute the FM polarization tensor, we begin by evaluating the following general outer product of fermionic 
Green functions: 


doj f d^k 


27r J (27r)3 


Go ( ia; + i^,k+ O Go | ) - Gq (iw,k) O Go (iw, k) 


= |q| 


<&!(§, q)i4 o I4 + ^ q)i4 o r„ 

n 

^ , q)r^ ^ ^ .") (p, q)r„ 0 F 


(Al) 


where the UV-divergent piece has been subtracted off, as described in the main text. Once this quantity is determined, 
any response function (and, in particular, the response function in the FM channel) can be calculated by contracting 
(Al) with appropriate matrices. The scaling functions appearing in (Al) are given by 


^ii§A) = 


d^x 

(27r)3 


(^x+q/2 4“ ^x—q/2) 


(ex-q /2 + ex+q/ 2 )^ - 


1 

4ex 


(A2) 








d?x 


2^2 


d„(x- 3) 


(27r)3 


^Stc — n 


q/2 


(£^x-q /2 + ^x+q/ 2 )^ 


(A3) 




d X J (^x+q/2 “1“ ^x—q/2)^m(^ “1“ 2^ 

(27r)^ 


2 ^x+q/ 2 ^x—q /2 


(£x-q/2 + ^x+q/2)^ ” 


4e3 


(A4) 


From (A4) it can be noted that < 1 )^'"’"^ = $ 4 ”’"*^ In the above equations we have introduced ci, 2,3 = ci and 04,5 = C 2 . 
In what follows it is useful to introduce matrices such that 


dm (^) — ^ ^ XiA.^ 
i,j = i 


(AS) 


where the matrices A"* are 


/ 0 0 0 
Al = — I 0 0 1 
2 ' 0 1 0 


A4 = 


/Q / 1 0 0 

' 0 ~1 0 

0 0 0 


/q / 0 0 1 

A" = ^ 0 0 0 

2 ' 1 0 0 


A3 = 


A4 / 0 1 0 

^ 0 0 
0 0 0 


A'l = 


-10 0 
0-10 
0 0 2 


(A 6 ) 


























These are proportional to the five symmetric Gell-Mann matrice^^ A^: 

2 


V3 


(A ,A ,A ,A ,A ) — (Ag, A 4 , Ai, A 3 , —Ag). 


(A7) 


We shall also make use of the fact that these matrices are traceless (^- A™ = 0). 

While it is possible to integrate the scaling functions (A2 )-(|A4|) numerically for a fixed value of the external 


momentum q, more convenient expressions can be obtained by separating out the q dependence from these integrals 
by applying the following transformation: 


Xi = -I- - qiqj) COST? -I- tijkqk smr]]x'j 

3 


(A8) 


and averaging the integrand over the angle 77 . Because this transformation rotates the integration variable x around 
the fixed unit vector q by an angle 77 , one has x'^ = x^ and q • x' = q • x. The following product of rotation matrices 
appears when (A8) is substituted into the scaling function integrals (A2)-(A4): 


^ ^ij^pq T i.^ip'^jq T TftpTjg) cos 77 4 “ {f^ipa^jqqa 4 ” ^jqbLipqt) Sin 77 

1 1, . . , . 

4 “ ^^ij^pq COS 277 4 “ ^\€-jqaT^ipqa 4 “ ejp^, 7 ^g 77 &) 

where we have introduced Lij = qiqj and T^- = 5ij — qiqj, and used 

^ ^ ^ipa^jqbqaqb — ^ij^pq ^iq^pj' 
ab 


(A9) 


(AlO) 


We have also used the fact that, due to the structure of the tensors and vectors with which L and T are contracted, 
we are permitted to interchange the indices i j, k I, p q, and r o s. (For this reason we use an arrow rather 
than an equality in (A9).) Contracting with the appropriate tensors and averaging over the angle rj, we obtain 


A^ix ± lUx ± = A-(x ± lUx ± I), ( ], 


(All) 


A™A^,(cr + f )p(a: + f ),(cr - f ),(a: - 


+ 2)p(* 2 ^q(^ ~ 2 ~ !)« ( ^ij^pq^klLrs 4" -LijLpgTklT, 


1 


-'ij kl rs 


4- -^LklLrsTijTpq + 2Lil~LprTjlTqs + -TijTpqTklTrs + -^TilTpsTjkTqr 


(A12) 


In obtaining this last expression, it was necessary to apply the angular rotation and averaging a second time in order 
to fully separate out the dependence on q. 


1. Case of isotropic fermion dispersion 

Beginning with the fully isotropic case (ci = C 2 and a = 0), using the rotation of the integration variable described 
above, and taking traces with the appropriate matrices, the scaling functions appearing in § are given by the 
following integrals over momentum: 


if) = (fp 


drr‘‘ < —TT — 


1 

4p2 


1 4- 4r^ 


(2r 


+ t) - ( xf 




II ifl / / . X 

l + 2^ + Jl-2^ > (A13) 


647r 
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dz 


{P" + + i) + \p^ + fp^ + h) _ + jp^ U 14 ) 


{P 


■2^Z2+ l)-_^2 


(2p>+222 + l)"-(f)' 


4(p2 -(- 2:2)3 


/4 ( ^ ) = 


(2^)2 7-c 


dz 


dp 


{p- + z^ + i) (2z^ - ip^ - i) 


(p2 + 22 + \y 


- 


(2p2 + 2z^ + l)"-(g)' 


2z^ - \p^ 
4(p2 + 22)3 


(A15) 


/5 ( S ) = 


(27r)2 7_, 


dz 


dp p^ 


{p^ + z^ + \) 


(p- 


2^z2 + i)"-z2 


(2p>+2,2 + 1)“-(f)' 


16(p2 + 22)3 


■.(A16) 


The integration coordinates are related to the original momenta k and q as z = k • q/q^ and p^ = k^/q^ — z^. 
In (A14)-(A16) it is possible to perform the integrals over z numerically, then it is straightforward to perform the 


remaining integrals numerically, as they are UV convergent due to the subtraction of 110 ( 0 , 0 ) in (|^. 

In the case of isotropic fermion dispersion (ci = C 2 ) but arbitrary Yukawa parameter a, the scaling functions 


determining the polarization tensor in (12 1 can be expressed as the following linear combinations of /i(^) in (AI3)- 


(A16): 


hi {^§, a ) =^(178/i + 22/3 + 22/4 + 267/5) - ^(190/i + SS/a + 68/4 + 286/5) cos(2a) 
+ i(82/i - 2/3 - 2/4 + 123/5) sin(2a). 


(A17) 


/i 2 (^,a) = ^( 26/3 + 17 / 4 ) [cos(2q!) - 1] - 9(/3 + /4)sin(2a), (A18) 

01 

d 3 (^,a) = ^/3 [1 - cos(2a)], (A19) 

= |t(404/3 + 743 / 4 ) - |-(92/3 + 173 / 4 ) cos(2a) + ^( 28/3 + 67 / 4 ) sin( 2 a), (A20) 

^64 64 8 

^ 5 (^, 0 ) = ^/a [1 - cos(2a)] + 5 / 3 sin( 2 a). (A21) 


2. Fully anisotropic case 


In the more general anisotropic case where ci 7 ^ C 2 , we introduce the anisotropy parameter 5 for the fermion 
dispersion by letting ci _2 = 1±(5. The integrands in (A2)-(A4) can then be expanded perturbatively to leading order 
in 5. The fermion dispersion becomes 


x±k/2 


= (x±i) 


2 , +E™=lCmC(x± I) 


(X±|) 


0(^2), 


(A22) 


wher e we h ave introduced = (!> 1,1, —1, —l)m- A s in the isotropic case, we can separate out the dependence on 
q in (A2)-(A4) by repeatedly applying the rotation (A 8 ) and averaging over the angle 77 . (In some cases products 
of more than four rotation matrices appear, leading to a much longer and more tedious calculation in which the 


transformation (A8) must be performed up to four times.) This procedure results in the following expressions: 


$1 +5 +$ 1 , 3 ( 1 ?) ^/a(q-A“.q )2 


0{6% 


{A23) 
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(f.«) 


=* 


2,1 


f )q-A".q + (5 


$2.2 (f) q • A" • q + $2,3 (f) q • A- • q^^a(q • A“ • q)^ 


+ $ 2.4 (f) E • A“ • q)(q • A“A" • q) + $2.5 (f) E ^ 


(A24) 


+ 0{S% 


$ 


(f , q) = $4.1 (f) (q • A’" • q)(q • A" • q) + $4.2 (f) q • A'"A" • q + $4.3 (f ) Tr[A™A"] 


+ (5 


<1 (§) + (fill (q ■ A™ • q)(q • A" • q) 


4,4 I q: 


+ 


<1 (§) + i^rn + e„)$i'l f q • A”^A" • q 


4,5 \ q2 


+$4.10 (i 


+ $4.12 ( 3 


+ $4.16 ( 3 


{§) + i^rn + en)$i!^ (f)] Tr[A™A"] + $4.7 (f) (4 ' A™ A" • q) ^^^(4 ' A“ • 4)" 

a 

+ ‘I’ 4.8 (f) (4 • A" • 4) E• A“ • 4)(4 • A“A- • 4) 

a 

+ $ 4.9 (f) (4 • A™ • 4) E • A“ • 4)(4 • A“A- • q) 

a 

(q • • q)(q • A“A" • q) + $4,n ^ ea (4 ' A“ • 4)(4 • A™A“A” • q) (A 25 ) 

a 

(4 • A“ • 4)(4 • A™A”A“ • 4) + $4.13 (f) 51 • A“ • 4)(4 • A^A^A- ■ 4) 

a 

+ $4.14 (§) (4 • A" • 4) 5] ?a(4 • A“A”^A“ • 4) + $4.15 (f) (4 • A™ • 4) 5] Ca(4 ' A“A"A“ • 4) 

a a 

q • A“A'"A”A“ • q + $4.17 (^) 51 ■ A™A“A”A“ • q 

A-.q)(q.A".q)5]Ca(4-A“-4)" 


+ $4.18 ( S 


+ ^4.20 ( 3 ] 


,q • A"A“A'"A“ • q + $4.19 
"A"] 5] Ca(4 • A“ • 4)" + $4.21 (f) 5] CaTr[A™A’^A“]q • A“ • q I + 0{S3. 


The scaling functions appearing in (A23)-(A25) are similar in form to those given in the isotropic case, though we 
shall not list them all here explicitly for the sake of brevity. 


Once the scaling functions are known, taking the traces in (141 gives the following expressions for the components 
of the polarization tensor: 


n“(zfi,q)-n5"(o,o) = 


32' 


(328sin(2a) - 285 cos(2a) + 445)$i(^, 4) 


+ (120sin(2a) + 21cos(2a) + 75) ($f’^^(p,q) + $^"’^4||,q) - 


s(3’3) / iQ 


s(i>i)nn 


+ (216sin(2a) - 219cos(2a) + 315)$f’'‘^(||,4) 

+ V3(-224sin(2a) + 132cos(2q;) - 260)$^^^^(f|, 4) 


(A26) 


+ (-8sin(2a) -87cos(2a) + 55)$f’®^(p,4) 


n7(7f],q)-n7(o,o) 



3(72sin(2Q;) - 73cos(2a) + 105)$i^’^^(p, 4) 


+ 2y3(-80sin(2a) + 75cos(2a) - 107)$f4) , 


(A27) 
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n-(zii,q)-n-(o,o) 



(216sin(2a) - 219cos(2a) + 315)$i^’^^(^, q) 

+ (240sin(2a) - 225cos(2a) +321)4>f’'‘^(p,q) 

+ 73(80sin(2a) - 75cos(2a) + 107)$f(p, q) , 


(A28) 


n^ 7 *f^,q)-n™(o,o) 



(328 sin(2Q;) 


285cos(2a)+445)$i(^,q) 


+ 3(40sin(2a) + 7cos(2a) + 25) "^(f ,q) + #"^(f, q)) 

+ (216sin(2a) - 219 cos(2a) + 315)$i‘‘’'‘^ ((p, q) (A29) 

+ 473(56sin(2a) -33cos(2a) + 65)$f’^P(p,q) 

+ (-8siii(2a) - 87cos(2a) + 55)4>fq) , 


n^/(*ii,q)-n^/(0,0) 



(-240sin(2a) + 225cos(2a) - 321)$^^’^^^, q) 


+ 73(80sin(2a) - 75 cos(2a) + 107)$^^’®^^, q) 
+ (216sin(2a) - 219cos(2a) + 315)4>fq) , 


(A30) 


n;^(*o,q)-n;^(o,o) = -|q 


(328 sin(2a) - 285 cos(2a) + 445) ($i (f, q) + , q)) 

+ (120 sin(2a) + 21 cos(2a) + 75) , q) + (^, q) 

-$rHi,q)-#"^(fl,q)' 


(A31) 


The remaining off-diagonal components follow from 11^ = 11;^ . 
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